Two-loop QCD corrections to top quark decay 
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We present a determination of a new class of Feynman diagrams relevant for second-order QCD 
corrections to the top quark decay t — > bW. Modern computing techniques allow us to perform a 
reduction of the original loop integrals to master integrals. We obtain the analytical value of the 
top decay rate as an expansion around the limit of massless b and W. 



I. INTRODUCTION 

It is well-known that the top quark is the heaviest 
fcrmion in the Standard Model. An immediate conse- 
quence of that fact is its very short lifetime which turns 
out to be shorter than the characteristic timescale of non- 
perturbative QCD effects. Thanks to this feature, the 
top behaves almost like a free quark, which is really a 
unique situation in QCD. Moreover, the process t — > bW 
is an absolutely dominant decay channel. Additional de- 
cay channels arise in various extensions of the present 
Standard Model (like t — + bH or t — ► supersymmetric 
particles). For these reasons, top physics seems to be a 
particularly clean laboratory for new physics searches but 
two complementary pieces of information are necessary. 
First, we need an accurate experimental determination of 
the top decay rate, which will be feasible in the future, 
e.g. at the Next Linear Collider. Also, to fully utilize this 
measurement it is incredibly important to have theoreti- 
cal predictions in the framework of the present Standard 
Model with accuracy matching the experimental preci- 
sion. 

Much effort has been invested into studies of radiative 
corrections to the top quark decay width. The one-loop 
QCD contribution was computed in 0. It decreases the 
tree-level rate by about 8.4%. The electroweak NLO part 
was evaluated in 2] but its effect turned out to be much 
smaller — it increases the decay rate by about 2%. More- 
over, the electroweak contribution is almost entirely can- 
celed when one takes into account the finite W width [3j . 
The NNLO QCD contribution has been studied so far 
in various kinematical limits. For example, there exists 
a result in zero recoil kinematics with m\y set to zero 
, as well as a numerical study done by means of Pade 
approximations ||. All these calculations estimate the 
NNLO QCD contribution at about —2% but, for compar- 
ison with planned measurements, a precise and preferably 
analytical result would be of interest. 

In the present study we start with the limit of a mas- 
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sive top quark decaying into a massless b quark and a 
massless W boson. Then we construct an expansion in 
the artificial parameter (mw /mt) 2 ■ Taking sufficiently 
many expansion terms, we can closely approach the phys- 
ical point (mw/mt) 2 ~ 0.213 and give a precise predic- 
tion for the second order QCD contribution to the decay 
rate. 



II. METHODS USED IN THE CALCULATION 

Let us now briefly discuss the computing methods 
which were used in the calculation. In the traditional 
approach the contribution to the decay rate can be di- 
vided into three distinct classes: real radiation with glu- 
ons in the final state, virtual loop corrections, as well 
as a mixed real- virtual case. Effective tools for comput- 
ing multi-loop diagrams have been developed over the 
last few years. Conversely, real radiation requires inte- 
gration over phase space in the final state, which has 
to be done manually on a case-by-case basis. This usu- 
ally constitutes the bottleneck of the computation and an 
automated method for dealing with real radiation seems 
mandatory at NNLO. 

A very simple but powerful idea is to map all real radi- 
ation diagrams onto loop integrals. To achieve this goal 
we apply the optical theorem, which relates contributions 
to the decay to imaginary parts of self-energy diagrams. 
In Fig. ^ we show three examples of the diagrams that we 
have to consider in order to calculate t — > bW at 0(a 2 s ). 
There are 19 abelian diagrams, 11 non-abelian diagrams, 
and 6 diagrams with quark vacuum polarization loops 
that have to be taken into account. The various cuts 
of these diagrams correspond to two-loop virtual correc- 
tions or emissions of one or two real quanta. However, 
the price to pay for the automation gained by expressing 
these diagrams in terms of a uniform set of self-energy 
graphs is the necessity of computing the imaginary parts 
of three-loop diagrams. This difficult task can be orga- 
nized in the following way. First, we perform a reduction 
of the original integrals to a set of simpler master inte- 
grals. In this phase we use the so-called integration-by- 
parts identities^ which lead to a system of recurrence 
relations between various loop integrals. The solution of 




FIG. 1: Examples of diagrams whose cuts contribute to the 
0(o? s ) decay rate for t — * bW: (a) light or heavy quarks; (b) 
abelian; (c) non-abelian. 



such a system provides reduction formulas which allow us 
to express the original loop integrals in terms of master 
integrals. 

There are two main paradigms for solving the system 
of recurrence relations. In the traditional method, one 
inspects the structure of the identities and rearranges 
them manually for an efficient iterative solution of the 
system. This "by inspection" method has proven to be 
very successful in numerous applications (see, e.g., - 
but requires much human effort to implement. A 
newer method is based on the Gaussian elimination of a 
large system of identities with a given orderin g fu nction 
which measures the "difficulty" of each integral 11]. This 
approach is fully automated and process-independent, 
but requires significant computing power. Thanks to an 
effective implementation and the growing performance 
of hardware, this method has become feasible at NNLO 
quite recently. In the present project we used both meth- 
ods in parallel. We obtained identical results which serve 
as a crosscheck for the calculation. It also allows us to 
compare the two approaches and point out strengths and 
weaknesses of both. 

The final step is the evaluation of the master integrals 
arising in the problem. For t — ► bW at 0{a 2 s ) there are 23 
master integrals which have been computed analytically. 
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X corresponds to the Born decay rate and X\ is the 
already known NLO correction. Our goal in the present 
calculation is the coefficient X 2 , which can be subdivided 
into four gauge-invariant color pieces: 

X 2 = C F (T R N L X L + T R N H X H + C F X A + CaX na ) , 

(3) 

where Cf = 4/3, Ca = 3, and Tr = 1/2 are the usual 
SU(3) color factors and Nl and Njj denote the number 
of light (m q = 0) and heavy (m q — m t ) quark flavors. 

We obtained series coefficients to uj 5 analytically. In 
terms of numerical values our expansion reads: 

X 2 = -12. 579+3. 743w+6.054w 2 +5.619w 2 logw+O(w 3 ). 

It is a straightforward task to improve the accuracy of 
an expansion by computing more terms in to. Moreover, 
the present expansion can be smoothly matched with the 
one in the zero recoil kinematics studied previously |9j in 
the context of semileptonic b quark decays. The result 
of such a matching procedure is depicted in Fig. [21 
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FIG. 2: Matching of expansions around u = (thick line) 
and (j = 1 (thin line). The solid line denotes the resulting 
decay width valid in the full range of uj. 



III. RESULTS 

Taking advantage of the tools discussed in the previous 
section, we can obtain the decay rate in the form of an 
expansion around the limit mw = Wty = in the param- 
eter 10 = (mw/mt) 2 (see 0). The decay width can be 



For the measured ratio of the W and top masses |13j , 
uj ~ 0.213, our expansion gives X 2 = —15.5(1), where 
the uncertainty is almost entirely due to the experimental 
error in the determination of m t . The theoretical error, 
which originates from taking a finite number of terms in 
our expansion, is 20 times smaller and can be reduced 
further if needed. Using a s (m, t ) — 0.11, we find that the 
two-loop correction decreases the tree level decay rate by 
about 2%, in agreement with earlier expectations. 
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IV. CONCLUSIONS 

We have presented a new analytic 0(a 2 s ) result for the 
decay t — » bW in terms of a parameter to = (mn//m t ) 2 
and in the limit of mj = 0. This result has enabled us 
to confirm or modify slightly the corresponding results of 
previous numerical calculations. Our formulas are readily 
applicable to other physical processes such as muon decay 
and the semileptonic b quark decay b — > ulv. 

Our results depend on the imaginary parts of a novel 
class of three-loop integrals, which we have obtained us- 
ing two independent paradigms for the solution of large 
systems of recurrence relations. To the best of our knowl- 



edge, this is the first time that both approaches have been 
used simultaneously to obtain a new result, and an ob- 
jective analysis of the strengths and weaknesses of each 
approach will increase the efficiency of other large calcu- 
lations in the future. 
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